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Abstract. We show that many features of the recent experiments of Schneble 
et al. [D. Schneble, Y. Torii, M. Boyd, E.W. Streed, D.E. Pritchard and W. 
Ketterle, Science 300, ^75 (2003)], which demonstrate two different regimes 
of light scattering by a Bose-Einstein condensate, can be described using a 
one-dimensional mean-field quantum CARL model, where optical amplification 
occurs simultaneously with the production of a periodic density modulation in 
the atomic medium. The two regimes of light scattering observed in these 
experiments, originally described as "Kapiza-Dirac scattering" and "Superradiant 
Rayleigh scattering", can be interpreted as the semiclassical and quantum limits 
respectively of CARL lasing. 



PACS numbers: 42.50.Fx 

1. Introduction 

The study of nonlinear optical phenomena arising from the collective motion of 
atoms in dynamic optical fields has been an active field of predominantly theoretical 
research over the last decade. A large fraction of this work has been concerned 
with the Collective Atomic Recoil Laser (CARL) Q] [0 El ED E] As part of the 
continuing progress in the production and investigation of ultracold atomic gases and 
Bose-Einstein condensates (BECs), there have been several experiments which have 
demonstrated the validity of these models and realized some of their predictions. 
Examples include the observation of collective atomic recoil lasing by a cold thermal 
gas in a high-finesse cavity |H] and the observation of superradiant Rayleigh 
scattering by a BEC [Q]. Superradiant Rayleigh scattering involves the 
production of pulses of coherently scattered radiation simultaneous with the splitting 
of the condensate into discrete momentum groups due to atomic recoil. Several 
theoretical models have been used to describe the evolution of the superradiant 
scattering process El D2I El UM including an extension of the original 
classical CARL model 00 0] to include a quantum treatment of the atomic dynamics 
[SJElEl- Recent experimental work by Schneble et al. has shown that in addition 
to the superradiant Rayleigh scattering process originally observed in 0, there is 
a second scattering regime termed "Kapiza-Dirac scattering". During superradiant 
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Rayleigh scattering, the scattering process involves only emission of scattered photons 
i.e. absorption of pump photons and emission of scattered (probe) photons. In 
contrast, during Kapiza-Dirac scattering the scattering process involves both emission 
and absorption of scattered (probe) photons i.e. absorption of probe photons and 
emission of pump photons. 

In this paper it is shown that many features of the recent experiments of 
Schneble et al. can be described using a one-dimensional mean-field quantum 

CARL model, where optical amplification occurs simultaneously with the production 
of a periodic density modulation in the atomic medium. Using this model, we 
demonstrate that the two regimes of "Kapiza-Dirac scattering" and "Superradiant 
Rayleigh scattering" observed in can be interpreted as the semiclassical and 

quantum limits respectively of CARL lasing. It will be shown that the two regimes 
are distinguished by the relative size of the gain of the superradiant scattering process 
and the frequency separation of the absorption and emission peaks. A significant 
difference between the results presented here and those of other theoretical models 
of the experiments of Schneble et al. E| ls that in this model the regime of 
scattering is not determined by the pump pulse duration, so the scattering process 
does not evolve in time from one to the other. 



2. Model 

The model used to describe the BEC-light interaction is the mean-field quantum- 
CARL model originally derived in |SJ El ■ The model is one-dimensional and describes 
the evolution of a backscattered (probe) field arising from scattering of a pump laser 
field (assumed to be of constant amplitude) by an elongated BEC. In the experiments 
of Schneble et. al. the geometry of the experiment is essentially two-dimensional, 
as illustrated schematically in fig. [3 with emission of two endfire modes from each end 
of the long axis of the condensate, propagating transversely to the pump laser. If we 
assume that coupling between the endfire modes (which are much weaker than the 
pump) is negligible, each endfire mode can be assumed to evolve independently and 
the atomic motion is one-dimensional. 

When the pump laser is sufficiently detuned from the atomic resonance, it leaves 
the atoms in the internal ground state. Consequently radiation pressure due to 
absorption and subsequent random incoherent, isotropic emission of a photon, can 
be neglected. In this detuned regime, coherent scattering of the pump laser is the 
dominant process. The atoms interact with a laser beam of wave vector k and 
scatter photons of wave vector k s , recoiling with a momentum fiq = %{k — k s ). The 
atoms, initially scattered randomly into various momentum states, interfere with the 
atoms in the original momentum state. This creates a matter wave grating having 
the correct periodicity to further scatter the laser beam in the direction k s . In an 
elongated condensate a preferential direction for the scattered photons emerges causing 
superradiant Rayleigh scattering. Both the matter wave grating and the scattered light 
are coherently amplified [HI E3 • 

In a simplified ID description of the process along the direction of the atomic recoil 
momentum Tiq, the evolution of the matter wave field \&(#, t) and of the dimensionless 
amplitude a(t) of the scattered radiation is determined by the following quantum 
mean-field CARL model 00: 

.5* 9 2 * 
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Figure 1. Schematic diagram illustrating the geometry of the superradiant 
scattering experiments of 1161 151 
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where 9 = qz (with q = \q\ w V%k), a = (eoV/2fiu} s ) 1 / 2 E is the dimensionless 
electric field amplitude of the scattered beam with frequency u> s , u> r = hq 2 /2m is 
the two-photon recoil frequency, g = (Q/2A)(ujd 2 /2h~€oV) 1 / 2 is the coupling constant, 
O = dEp/H is the Rabi frequency of the laser field with constant amplitude E p and 
frequency to — ck, A = uj — loq is the detuning from atomic resonance ujq, d = e • d is 
the electric dipole moment of the atom along the polarization direction e*of the laser, 
V = AL is the volume of the condensate, A is its cross-sectional area, L is its length, 
N is the total number of atoms, and S = to — oj s . The matter wave field is normalized 
such that f_°° d9\^li\ 2 = 1. The second term on the right hand side of Eq. JJJ is 
the self-consistent optical lattice, resulting from the interference between the laser 
and the scattered radiation, whose amplitude is amplified by the matter wave grating 
described by the first term on the right hand side of Eq. (J2J. The third term on the 
right hand side of Eq. JTJ) describes the mean-field effect of the atom-atom interaction 
due to binary collisions, where ft = 8nh~qa s N/mA and a s is the scattering length. Eq. 
P]l has been written in the 'mean-field' limit, which models propagation of light with 
respect to the atoms by replacing the non uniform amplitude by its average value and 
by adding to the equation a damping term with a decay rate k ~ c/2L of the order of 
the inverse of the photon flight time through the condensate. 

Notice that the interaction time t appearing in Eqs. JJJ and (J2J is the pump 
field duration, so that our model is suitable to explore both Superradiant Rayleigh 
scattering ('long pump pulse duration') and self-stimulated Kapitza-Dirac diffraction 
regime ('short pump pulse duration') investigated in ref. |16j . It will be shown here 
that the relevant parameter in these experiments is not the pump pulse duration but 
the characteristic time of the optical/matter- wave amplification process, given by the 
inverse of the gain rate. 

If we assume that the atomic wave function is periodic on the scale of the optical 
potential, with spatial period 2n/q = A/V2 which corresponds to a period in 9 of 2n, 
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then we can expand the atomic wave function in a Fourier series 

oo 

*(0,t)= cn(t)e m ( e +^ (3) 

n— — oo 

Furthermore, we assume that the condensate is sufficiently dilute such that 
4:Trha s n s /m -C u> r , where n s is the average atomic density, so that the atom-atom 
interaction term in Eq. may be neglected. The effect of the atom-atom term on 
the collective recoil lasing has been investigated in ref. |18j . 

Substituting for 9(6, t) using Eq. ©, it can be shown [SJEIE] that Eqs. (QJ and 
J2J can be rewritten as 
dc n 

= - in(uj r n + S)c n - g (ac„_i - a*c n+1 ) (4) 
d °° 

^ ^ c » c »-i - Ka > ( 5 ) 

n— — oo 

In this quantum description, the Fourier expansion is equivalent to expanding the 
wave function 9(9,t) in the set of momentum eigenstates with eigenvalues p — (hq)n 
and p n = \c n \ 2 is the probability for an atom to have a momentum p — (Tiq)n. 

In the superradiant regime explored in the experiments of ref . [9"1 ITUl 111! I16| . the 
radiation damping rate n is always much larger than the gain rate and/or the recoil 
frequency u) r , so that the field amplitude follows the atomic motion adiabatically . 
Hence, neglecting the time derivative, Eq.JSJ) yields 

a re > c n c n _ x . (6) 

K < 4 



3. The semiclassical and quantum limits of the superradiant regime 

In order to obtain the gain coefficient of the superradiant process in the semiclassical 
and quantum limits from the dynamical equations, let us consider the initial 
equilibrium state with no field, a — 0, and all the atoms at rest, i.e. in the 
momentum state n = 0, with cq = 1 and c m = for all m 7^ 0. Linearizing 
around this equilibrium solution and neglecting the small detuning S between the 
laser and scattered frequencies, Eqs. (@J and I© reduce to the single linear equation 
for B = c 1 +c*_ 1 : 

d 2 B 

— Y +u r (u r -iG)B = 0, (7) 
where G is the superradiant gain 

G = ^lE = ^ TT scN 

r is the natural lincwidth, r sc = r(17/2A) 2 is the Rayleigh scattering rate and 
I sat — ceh 2 T 2 1 'Ad 2 is the saturation intensity. Eq. © coincides with the superradiant 
gain coefficient reported in Eq.(4) of 

It is easy to show that Eq.J7|) has an unstable solution B(t) oc exp[(Ai + iAa)t], 

with 

1/2 



■Jr^Ju) 2 + G 2 
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Ai = ' Vr r 0) 
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(10) 



|c-i| y/iXx + G/2) 2 + (A 2 + u r y 



These expressions show that the instability may have a semiclassical or quantum 
character, depending on the ratio between the superradiant gain G and the recoil 
frequency uj t . In fact, in the limit where G <C ui r) Eqs. 10 and 1|10[) give Ai « G/2, 
A2 ~ u) r and |ci|/|c_i| « G/(4ov) <C 1. In this limit, the SR process is quantum in 
nature, with only the lower state n = — 1 being populated, and the SR gain is G [5lll2|. 
In the opposite limit in which G » ui r , Eqs. © and (fTTTf) give Ai k A 2 ~ ^Juj r G/2 and 
l c i|/l c -i| ~ 1 — \/2u) r /G. In this limit, the SR process is semiclassical in nature, with 
the states n = 1 and n = — 1 almost equally populated. In this case, the superradiant 
gain is given by 0] 



The SR gain in the semiclassical limit is always lower than the SR gain in the quantum 
limit, since G'/G — y/2uj r /G <C 1. Notice also the different dependence of the SR 
gain on N in the quantum and the classical limits. 

4. Numerical Results 

In order to observe the behavior of the collective scattering process in the nonlinear 
regime, Eqs. (0J and JSJ) were integrated numerically with initial conditions a = 0, 
c_i = l/y/N, Co = \J\ — l/N, and c m — when m 7^ —1,0 and with parameters 
corresponding to the semiclassical and quantum regimes of evolution. These 
parameters correspond to the experiments of ^Sj i-e. a 87 Rb condensate illuminated by 
a pump beam of wavelength A = 780nm and an intensity of 63mW cm -2 . The pump 
couples to the 5Si/ 2 — > 5P 3 / 2 transition which has a natural width T — 0.37 x 10 8 
s _1 , dipole moment d = 2.07 x 10~ 29 C-m, saturation intensity I sa t = 2.5mW cm -2 
and recoil frequency uo r = 4.7 x 10 4 s _1 . The condensate had a cigar-shaped form, 
15 /xm in diameter and 200/im in length, so that k = 7.5 x 10 11 s" 1 . Using these 
parameters, the quantum and classical superradiant gain coefficients of Eqs. (jHJ and 
JllJ are G w 4.9 x 10 6 x 7V/|A| 2 and G' w 6.8 x 10 5 x VN/\A\, respectively, where A 
is the pump-atom detuning in MHz. We assume that S — and that half of the atoms 
in the condensate participate in each of two superradiant emissions along the main 
axis of the condensate and that the number of atoms participating in the collective 
scattering process is N = 10 5 rather than N = 10 6 as quoted for the number of atoms 
in the condensate in Qualitative support for this assumption is provided by fig. 3 A 
of m3 , which shows a large fraction of the condensate atoms do not participate in the 
coherent superradiant scattering process. From |16| . the 'Kapiza-Dirac' experiment 
was carried out using A = — 420MHz, so that g = 3.2 x 10 6 s _1 and the quantum 
superradiant gain is G = 3 x 10 6 s _1 . Consequently the collective scattering process 
is semiclassical in nature, since G/uj r ~ 58. For the experiments in the 'superradiant 
Rayleigh scattering' regime, A = —4400MHz, g = 3.07 x 10 5 s _1 and the quantum 
superradiant gain coefficient is G = 2.5 x 10 4 s _1 , so that the collective scattering 
process is quantum-mechanical in nature, since G/oj r ~ 0.53. 

Fig. shows snapshots of the atomic momentum distribution at different times 
for the semiclassical case with A = — 420MHz. It can be seen from fig. EI a) that, in 




(11) 
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agreement with our theory, at the beginning of the interaction momentum states n = 1 
and n = — 1 are the only non-zero momentum states to have significant population, 
with pi ~ p_x, where p n = |c„| 2 . Eq. (jTU|l predicts a ratio |pi|/|p_i| ~ 0.7. As time 
progresses, many momentum states with both positive and negative n are populated. 
However, the momentum distribution is not symmetric about n — 0. The average 
momentum is less than zero, so from momentum conservation there is a net gain of 
the scattered radiation field, shown in figEl due to the difference in photon absorption 
and emission rates. Fig. shows that amplification of the scattered field in the 
semiclassical case is simultaneous with the growth of a strong density modulation in the 
condensate, as represented by the bunching parameter, b = (e~ l8 ) = c n c* l _ 1 . The 
atomic momentum distributions shown in fig.|2]are similar to the time of-flight images 
observed for the so-called "Kapiza-Dirac scattering" observed in ^] , where population 
of momentum states due to absorption and emission of radiation was observed. 
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Figure 2. Atomic momentum distribution in the semiclassical regime with 
G = 58oj r , when (a) t = llfis, (b) t = 15/is and (c) t = 17^ts. 



Fig. 01 shows snapshots of the atomic momentum distribution at different times 
for the quantum case with A = —4400MHz. It can be seen from fig. 0{a) that, 
in agreement with our model, at the beginning of the interaction momentum state 
n = — 1 is the only non-zero momentum state to have significant population, and that 
the population pi « 0. Eq. (|f Ofl predicts a ratio ~ 0.01. As time progresses, 

the atomic population moves sequentially from n = — 1 ^ n = — 2 — > ... and states 
with n > are never populated. This sequential decrease in atomic momentum gives 
rise to amplification of the scattered radiation field and again occurs simultaneously 
with the development of a strong density modulation in the condensate, as shown in 
fig. |S] In contrast to the semiclassical case, however, this amplification of the scattered 
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Figure 3. Flux of scattered photons, 2k|<j| 2 , and bunching factor, |b|, as a 
function of time in the semiclassical regime with G = 58aj r . 



field occurs due to emission of scattered (probe) photons only, and the spread in atomic 
momenta is much smaller than in the semiclassical case. The atomic momentum 
distributions shown in fig. 0] are similar to the time of-flight images observed for 
"Superradiant Rayleigh scattering" in where the atoms attain momentum only 

in the direction of the pump beam in discrete units of Tiq. 

5. Interpretation 

In |16| . the Kapiza-Dirac scattering regime and the Superradiant Rayleigh scattering 
regime are described as the "short-pump pulse" and "long-pump pulse" limits 
respectively. A problem with this classification is that it implies that for sufficiently 
long pump pulses the collective scattering makes a transition from Kapiza-Dirac 
scattering to superradiant Rayleigh scattering, or in the terminology of this paper 
a transition from the semiclassical to the quantum regime of CARL lasing. The 
results presented here however do not support this interpretation. They suggest that 
the distinguishing feature between the two experiments in JS] is not the pump pulse 
duration, but the atom-field detuning (which differs by an order of magnitude in 
the experiments of which determines the atom-field coupling and consequently 

the timescale of the superradiant scattering process. Rather than the ratio of the 
pump pulse duration relative to the two-photon recoil time (the inverse of the two- 
photon recoil frequency), it is the size of the superradiant decay time relative to the 
two-photon recoil time which determines whether or not the scattering consists of a 
sequence of emission processes, as observed in the "superradiant Rayleigh scattering" 
regime or simultaneous emission and absorption processes, as observed in the 

"Kapiza-Dirac regime" . 

A simple justification for this argument is as follows: The atoms are assumed 
to all be in the condensate initially, so that po = 1. The frequencies at which the 
probability of absorption [n — — > n — 1) and emission (n = ^ ra = — 1) 
are maximum are non-degenerate and differ by a frequency of 2cu r . A transition 
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Figure 4. Atomic momentum distribution in the quantum regime with G 
0.53u;r when (a) t = 0.3ms, (b) t = 1.0ms and (c) t = 1.7ms. 




Figure 5. Flux of scattered photons, 2k|<j| 2 , and bunching factor, |b|, as a 
function of scaled time in the quantum regime with G = 0.53oj r . 



from (n = — * n = — 1) will initiate superradiant or supcrfluorescent decay. The 
characteristic time of the superradiant decay is r sr ~ G _1 for k S> gvN j^J, so the 
spectral width of the superradiant pulse is ~ G — 2g 2 N/ k. If the spectral width of the 
SR pulse is much less than the absorption-emission frequency shift i.e. g 2 N/n < to r , 
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then absorptive transitions will not occur because they are non-resonant and only 
(n = — > n = — 1) transitions, i.e. absorption of pump photons and emission of probe 
photons, will occur. Consequently the system evolves in the quantum CARL limit 
when G < ui r . In contrast, if the spectral width of the SR pulse is sufficiently large 
(i.e. SR decay is sufficiently rapid) that it is much larger than the absorption-emission 
frequency difference, so that G 3> u> r , then absorptive transitions (n = — > n = 1) 
will be resonant and both emission and re-absorption of probe photons will occur. 
Consequently the system evolves in the semiclassical CARL limit when G ^> ui r - It 
should be noted that in contrast to previous explanations |161 H7) of the experimental 
results in 16 , in our argument the duration of the pump pulse is not a significant 
factor. The reason that the semiclassical CARL or "Kapiza-Dirac" regime can be 
observed using a short pump pulse in |16| is because the timescale of semiclassical 
superradiant decay, r' sr ~ 1/G" (see Ea. lfTT|) N | is much shorter than in the "Superradiant 
Rayleigh Scattering" example as a result of the decreased pump-atom detuning. In 
Ref. a suppression of quantum SR gain G of Eq. JSJ by around two orders of 
magnitude was observed and attributed in ^B] to the short duration of the pump pulse. 
The CARL model described here explains the reduced gain observed in JHj as a result 
of the fact that the SR scattering process is evolving semiclassically. Consequently 
the SR gain is not given by G as given by Eq. (jHJ but the semiclassical SR gain G' as 
given by Eq. (|1 1|) . The semiclassical gain, G' , gives a value consistently smaller than 
the one predicted by the quantum superradiant gain, G. Using the same parameters 
as those used in figs 0] and the ratio is G'/G » 0.16. 

6. Conclusion 

It has been shown that many features of the recent experiments of Schneble et al. > 
which show two different regimes of light scattering by the BEC, can be described 
using a one-dimensional mean-field quantum CARL model. The two regimes of light 
scattering observed in ^H] , described as "Kapiza-Dirac scattering" and "Superradiant 
Rayleigh scattering" in ^5], can be interpreted as the semiclassical and quantum 
limits respectively of CARL lasing. In the semiclassical limit, when g 2 N/ k ^S> uj r , the 
collective scattering process involves both absorption and emission of probe photons, 
and for sufficiently long times many momentum states are populated simultaneously. 
In the quantum limit, however, when g 2 N/n < u> r , the collective scattering process 
involves emission of probe photons only, and a maximum of two momentum states 
are populated at any time. We provide a simple explanation of these results in terms 
of a comparison between the frequency separation of probe emission and absorption 
events and the spectral width of superradiant decay between the initial and recoiling 
condensates. In contrast to previous models of the experiments in JHjj the pump 
pulse duration is not a significant factor in our interpretation. The results presented 
here support the view of the BEC-light interaction as a strongly coupled atom-optical 
system where the atoms and light (and consequently the matter-wave and optical 
gratings) evolve dynamically and self-consistently. This picture gives a more correct 
and complete description of the interaction than models in which the dynamics of 
either the atoms (matter-wave grating) or the optical fields (optical grating) are 
neglected. 
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